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It is anticipated that hard double parton scattering (DPS) will occur frequently in the collisions 
of the LHC, producing interesting signals and significant backgrounds to certain single scatter- 
ing processes. In order to make theoretical predictions of double scattering rates and properties, 
the double parton distributions (dPDFs) D-'p^^{xi,X2;Qa,Qb) are required. We discuss the first 
publicly available set of equal-scale (Qa — Qb) leading order dPDFs - the GS09 dPDFs. A brief 
account is given describing how pQCD evolution effects and sum rule constraints (the latter de- 
rived by us) have been incorporated into this set of dPDFs. We then present a summary of a 
phenomenological investigation into same-sign W pair production conducted using GS09. In 
this, the DPS signal produced using GS09 is compared with that obtained using simple products 
of single PDFs x (1 — jci — X2)", and the single scattering backgrounds {W^W^jj, di-boson and 
heavy flavour) are carefully calculated. It is found that the correlations in GS09 manifest them- 
selves in non-trivial kinematic correlations between the W bosons. However, it is unlikely that 
these correlations will be measurable at the LHC in the near future owing to the fact that the 
background is significant even after cuts. 
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1. Introduction 



Double parton scattering (DPS) is said to occur 
when two pairs of partons participate in hard interac- 
tions during a single proton-proton collision. This pro- 
cess is drawn schematically in figure [l|, in which all of 
the variables relevant to DPS have also been defined. 
Assuming the factorisation of the hard subprocesses A 
and B, and further assuming that the generalised dou- 
ble parton distributions ^ij{x\,X2,b;QA,QB) may be 
decomposed into a longitudinal piece and a piece de- 
pending only on the transverse separation b, we may 
write the cross section for DPS as follows: 
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(1.1) 

The function Dp (xi ,X2', Qa,Qb) is known as a double parton distribution function (dPDF) and 
may be interpreted at LO as the probability of finding a pair of partons in the proton with flavours / 
and j, longitudinal momentum fractions x\ and X2, when probed at scales Qa and Qb- The quantity 
Oeff is a non-perturbative factor which is related to the pair density distribution in transverse space. 
We take it to be a constant and equal to the CDF measured value of 14.5mb, with the caveat that the 
value of this quantity should be re-established at LHC scales using benchmark processes, m is a 
symmetry factor that equals 1 if A = B and 2 otherwise, and the are parton-level cross sections. 

In many of the extant studies of DPS (see [Q] for a complete list of references), a further as- 
sumption which has been applied is to take the dPDFs as products of sPDFs (i.e. ignore longitudinal 
momentum correlations). It is argued that this is justifiable at low due to the large population of 
partons at these x values, and there is some evidence from CDF ^ that it holds for sea partons at 
low X. Under this additional approximation, the DPS cross section reduces to a product of single 
parton scattering (SPS) cross sections: 



^DPS 



(1.2) 



We recall that the SPS cross section for a given hard process grows with energy as the single 



PDFs (sPDFs) are probed at lower x values where they are larger Equation ( |1.2| ) demonstrates that 
the DPS cross section goes approximately as the product of SPS cross sections - therefore DPS 
cross sections grow faster with energy than SPS ones, and DPS will be more important at the LHC 



than at any previous collider. It has already been established using ( |1.2| ) that DPS processes consti- 
tute important backgrounds to Higgs and other interesting signals (see e.g. In addition to this, 
DPS can be considered as an interesting signal process in its own right, as it reveals information 
about the correlations between pairs of partons in the proton. 

For these reasons, we need to ensure we have a good theoretical handle on DPS. It has already 
been established that the form ([T^ is not adequate. In [||] a LO 'double DGLAP' (dDGLAP) 
equation was derived describing the scaling violations of the dPDFs with a common factorisation 
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scale {Qa = Qb = Q)- A crucial prediction of this is that even if the dPDFs may be taken to be equal 
to products of sPDFs at one scale, then at any other scale the dPDFs will deviate from factorised 
forms. 

Further to this, simple factorised forms do not obey the following sum rule constraints, which 
we have shown are preserved by LO dDGLAP evolution if they hold at some starting scale [Q] : 

/ dxixxD^'^^{xi,x2;Q) = {l-x2)D^Hx2;Q) (1.3) 
h ^« 

^ ^juPp {^2\ Q) when 72 / ji or Ji 



Number Sum Rule: 



\-X2 



dxiDj'^j'-{xy,x2;Q)= < 



{Nj,^-\)D'p\x2;Q) when 72 = 71 
^(A^,-„ + l)D^^(x2;G) when72=7i (1-4) 



The symbol 71,, = 71 — (71 7^ g), and Nj^^, is the number of 'valence' 71 quarks in the proton. 
These sum rules have simple interpretations in terms of conditional probabilities. The first states 
that if you observe a parton with momentum fraction X2 in the proton, the momentum fractions of 
all other partons must add up to 1 — X2. The second states that if you observe a parton with flavour 
7 in the proton, the number of partons of flavour 7 elsewhere in the proton must be reduced by one 
(we use the term 'number effects' to describe this simple phenomenon). 

In this contribution we discuss the first publicly available set of LO equal-scale dPDFs - the 
GS09 dPDFs - which incorporates pQCD evolution effects plus the sum rule constraints. Section 
2 describes how the dPDFs were constructed incorporating these features [Q], whilst section 3 
summarises a phenomenological study of DPS performed using GS09 [S]. 



2. The GS09 dPDFs 

The GS09 dPDF package comprises a grid of dPDF values spanning the ranges 10-*^ <.xi < 1, 
10^^ < a;2 < 1, 1 GeV^ < < 10^ GeV^, which is available along with interpolation code from 
HepForge [^]. It has been obtained by constructing inputs that approximately satisfy the sum rules 
at 2o = 1 GeV, and then numerically evolving these inputs up to higher scales according to the 
dDGLAP equation. The sPDF set to which we have chosen our dPDF set to correspond is (a 
slightly modified version of) the MSTW2008LO set [^. 

Given the paucity of experimental data regarding the dPDFs, and in accordance with sim- 
ple arguments and the CDF results, we base our inputs on factorised products of MSTW2008LO 
sPDFs. However, we modify these basic forms in several ways to ensure that the input dPDFs 
approximately satisfy the sum rules. 

First, all of the dPDFs are multiplied by a factor p'j{xi,X2) which is designed to take account 
of phase space effects. This factor should ensure the appropriate behaviour of the dPDFs near the 
kinematic boundary xi + X2 = I - namely, a smooth decrease to zero. It was discovered that the 
following form for p'j gives rise to inputs which satisfy the momentum sum rules (plus appropriate 
number sum rules) well: 

p'j{xi,X2) = (1 -XI -X2)Hl -X2)-2-«« (2.1) 

a{i) is if / is a sea parton, and 0.5 if it is a valence parton. 
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We recall that there are only a finite number of valence quarks in the proton, as opposed to 
an infinite number of sea quarks and gluons. Number effects are therefore most significant in 
the context of valence quarks, and on this basis we have chosen to only take account of valence 
number effects in our inputs. This is done by dividing the UyUy part of any dPDF by two, and 
completely subtracting the dydv part. The reasoning behind this is that removing one up valence 
quark essentially halves the probability to find another, whilst there is no chance of finding two 
valence down quarks in the proton. 

Finally, we have added extra terms to input distributions whose flavour indices contain jj 
combinations to take account of so-called 'jj correlations'. These are essentially related to sea 
parton number effects, although they can alternatively be thought of as arising during evolution 
from some lower scale to Qq via g — jj splittings. It is important to include these terms in the 
equal flavour valence-valence (j'vj'v) inputs since the ' jj correlation' term is much larger for these 
than the quasi-factorised piece at low x. 

With these adjustments, our dPDF inputs satisfy all sum rules to better than 25% accuracy for 
a; < 0.8 (in the normal 'double human' basis). 

3. Effects of using GS09 dPDFs on same-sign WW DPS signal 

It is interesting to ask how the inclusion of proper pQCD evolution and sum rule constraints 
in GS09 affects experimentally measurable DPS signals. In [||], the same-sign WW DPS signal 
produced using GS09 was compared with that arising from simple factorised forms. The factorised 
forms used were simple products of MSTW2008LO dPDFs multiplied by (1 —xi —x-if, « = 0, 1,2 
(the 'MSTW,,' dPDFs). Same-sign WW production was chosen as the DPS process because it 
has been traditionally considered as a clean channel for observation of DPS. The cross section for 
same-sign WW production via SPS is suppressed to the same order of magnitude as the DPS cross 
section due to the large number of vertices required in the Feynman diagrams. What is more, this 
SPS background must always produce two jets in addition to the WW pair - so it can be efficiently 
removed via a jet veto. 

We chose to study the signal in the case when both Ws decay leptonically. A difference in the 
GS09 and MSTW,, sets can be observed in the following quantity which is particularly sensitive to 
longitudinal correlations: 



This quantity is plotted as a function of the minimum pseudorapidity cut on the detector hemi- 
spheres in figure The GS09 values for a^, are larger because the probability of a proton providing 
two large x (valence) partons is reduced under GS09 (due to number effects), so GS09 results in 
fewer events with the two leptons in the same hemisphere. 

In addition to a comparison of GS09 and factorised forms, [||] examines sources of background 
aside from the canonical SPS WW processes. It is found that there are several processes that can 
give rise to the same-sign lepton signal in the detectors: 

Heavy flavour: gg^tt, t^W^b^l^vb, t^W^b^qql^vc (3.2) 



g(?7/i X T];, < 0) - CJ{r||^ X T]/, > 0) 

a(T]/i X T7/2 < 0) + a(T]/i x tj/^ > 0) 



(3.1) 



gg^bb^BB + B^I+vX 



(3.3) 
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(a) Positively charged leptons 
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(b) Negatively charged leptons 



Figure 2: Pseudorapidity asymmetry uri, for pp collisions at = 14 TeV evaluated using different dPDFs. 
No cuts are applied. 



Electroweak gauge boson pair: 



qq' ^W+Z(f ) ^/+v/+(/") 

qq ^Z(f)Z(f)^ /+(/-)/+(/-) 



(3.4) 
(3.5) 



The leptons in brackets must fall outside the detector acceptance to give rise to the appropriate 
signal, which in our study meant that they had to satisfy > 2.5. The variety and strength of 
the backgrounds means that a careful choice of cuts is required to enhance the signal/background 
as much as possible - these are listed in [§]. Even with these cuts, the backgrounds are non- 
negligible and it is unlikely that a conclusive discrimination between GS09 and factorised forms 
will be possible using this channel in the near future. On the other hand, cleaner channels with 
regards to DPS may exist - we are currently examining double Drell-Yan as a possibility. 



References 

[1] CDF Collaboration, F. Abe et. ai. Double parton scattering in pp collisions at ^/s = l.8TeV, Phys. 
Rev. D56 (1997) 3811-3832. 

[2] A. Del Fabbro and D. Treleani, A double parton scattering background to Higgs boson production at 



the LHC, Phys. Rev. D61 (2000) 077502, [hep-ph/9911358] 



[3] R. Kirschner, Generalized Lipatov-Altarelli-Parisi equations and jet calculus rules, Phys. Lett. B84 
(1979) 266; V. P. Shelest, A. M. Snigirev, and G. M. Zinovjev, Gazing into the multiparton distribution 
equations in QCD, Phys. Lett. B113 (1982) 325. 

[4] J. R. Gaunt and W. J. Stirling, Double Parton Distributions Incorporating Perturbative QCD Evolution 



and Momentum and QuarkNumber Sum Rules, JHEP 03 (2010) 005, [arXiv: 0910 . 4347], 



[5] J. R. Gaunt, C.-H. Kom, A. Kulesza, and W. J. Stirling, Same-sign W pair production as a probe of 



double parton scattering at the LHC, arXiv:1003.3953 
[6] http://projects.hepforge.org/gsdpdf. 

[7] A. D. Martin, W. J. Stirling, R. S. Thorne, and G. Watt, Parton distributions for the LHC, 



arXiv: 0901 . 0002 



5 



